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ABSTRACT 

An approximate vacuum wave functioneil is proposed for 2+1- 
dimensional Yang-Mills theories. Using ^o. one can compute the 0"^""" 
glueball mass Mq in terms of the string tension. By using the idea of 
dimensional reduction, a prediction for Mq can be made in 3+1 
dimensions. One finds Mq =1.5 GeV. 
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I. Introduction 

Y£ing-Mills theories without quarks are expected to produce one or 
more bound states or glueballs. When qucirks cire present, such pure glue 
states should mix with qq states but some meson might have a dominant 
gluon content. Possible candidates for such states are fo(1300), fo(1590), 
fj(1710) and ti(1440).[1] The r|'(958) is expected to obtain a large fraction 
of its mass due to fluctuations in gluonic topological charge. [2-4] 

Initial lattice studies suggested that the lightest O"*""*" glueball might 
have a mass Mq of ~ 1 GeV. Such lattice gluebcill computations have been 
among the most numerlccilly demeinding due to a poor signal to noise ratio. 
However, due to adv£inces in gilgorithms and powerful computers, much 
progress has been made. The most recent lattice results give v£ilues of Mg of 
1550 ± 50 MeV [5] and 1740 ± 71 MeV [6], suggesting that the fj(1710) or 
fo(1590) might be the lightest O"""" glueball state. The tj(1710) is favored 
since it has decay widths consistent with numericcil simulations. [7] 

The purpose of this letter is to obtain an approximate analytic 
computation of the 0"^"^ glueball mass in terms of the string tension o. The 
calculation is carried out in D=2+l dimensions and extended by dimensionsil 
reduction to D=3+l dimensions. Anal5^ic studies often provide more 
physical insight than numerical simulations. In addition, the interplay 
between numerical and analytic approaches Ccin assist either approach in 
obtaining new methods cind results. 

Our method consists in postulating the form of the ground state wave 
functioneil ^o- The proposed is adjusted to agree with wesik £ind strong 
coupling limits. Support for our comes from the previous work [8,9]. In 
pcirticulcir, numerlccil cind aneilj^ic studies [10-17] of the lattice Hamiltonian 
formulation [18] give rise to vacuum functionals simileir to our ^o- 

The Hamiltonian for the D=d+ 1 dimensional Yang-Mills theory is 




(1.1) 
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where Ef(x) = Ig^^ - Bfj = d^Af - BjAf + f^""^^, g is the gauge coupling 

constcint cind f^^^ cire the structure constcints: [A,^,A,^] = if^^^A,^, where the Lie- 
algebra generators are normalized so that Tr{A^) = 5^^/2. 

In 2+1 dimensions, our key result is 



Mg-,,27v . (1-2) 



8c 
g2Cf ' 

where Cf is the v£ilue of casimir operator for the fundcimental representation: 
= Cf I. For SU(N], Cf = (N2-1)/(2N). 



n. Hie ./^proximate Vacuum Wave Functional 

Express the ground state as 

^0 = exp(-f(B)) , (2.1) 
where f is a functional of B. As g goes to zero, the perturbative fpt(B) is 

fpt(B) =^jd''xB;(x)[^-^B;(x)j 

-1/2 

where the kernel (-82) (x,y) is 



^2 
J 



IV , _ f d'p exp(ip • (x-y)) 



In the abelian case for which there are no interactions, perturbation theory 
cind Eq.(2.2) £ire exact. For the non-abelicin case, Eq.(2.2) is not gauge- 
inv£iriant but the violations of gauge invariance are of order g2. 

It has been conjectured that, in strong coupling, the ground state wave 
functional is governed by fsc(B) with [9] 

fsc(B) = |jddxB^Bj.(x) , (2.4) 
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where |i is a pcirameter with dimensions of length There is much 
evidence in support of Eq.(2.4).[9-17] In particulcir, lattice theory gives 
Eq.(2.4) as the leading strong coupling result. [15] 

In a confining theory in space-time dimensions D with 2 < D < 4, 
Eq.(2.4) can be derived by the following argument. Work in a basis of 
generalized closed Wilson loops of arbitrary shape and number and in 
arbitrary representations. This gauge -invariant set forms a complete set of 
variables. Expand the vacuum functional in such a basis. Using constructive 
field theory, the vacuum state can be obtained by doing a functional integral 
weighted by exp(-S) in which one integrates over half of space-time 
corresponding to t<0 and uses free boundary conditions at t=0. The 
functional at t=0, obtained in this way, is the exact ground state wave 
functional. In a confining theory, Wilson loops of large area arc suppressed. 
Due to the kinetic energy term g^jd^x E^ E^(x)/2 in H, one sees that the 

there is cilso a contribution per unit length to the energy. Hence, Wilson 
loops of Icirge area or large perimeter are suppressed. One concludes that, 
in a confining theory, f(B) in Eq.(2.1) is a sum over girbitrsiry numbers of 
Wilson loops in £irbitr£iiy representations for loops which are of smeill size. 
As a consequence, f(B) acts like a loccilized field theory. In computing 
vacuum expectation values, a confining theory in D dimensions becomes a 
localized field theory in d=D-l dimensions governed by an action equal to 
2f(B) (the factor of two cirises because l^o'^ enters in computations). Now, 
use the idea behind the renormsilization group. To compute the behavior of 
a Icirge Wilson loop or long-distcince correlation function, one may integrate 
out short-distcince degrees of freedom. Integrate out to a scale slightly 
larger that the confinement distance. Then, the resulting field theory will 
be dominated by the local gauge-inv£iri£int operator of lowest dimension. 
This operator is By By. In space-time dimensions D with D < 4, strong 

coupling corresponds to Icirge distcinces. Hence, the effective vacuum 
functional in the strong coupling limit is given by Eq.(2.4). 

From the above discussion, it is clear that the true vacuum functional 
is quite complicated. However, a simplified functional, which interpolates 
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between the wccik coupling (Eq.(2.2)) and strong coupling (Eq.(2.4)) forms, 
might produce good results for computations. Consider the approximate 
vacuum functioneil governed by 



-1/2 



«B) = ^\d'xTr[B^^iK){-'D^^^ + ml) (x,y)B^j(y) 



(2.5) 



where By(x) = A,^By(x), Tr stands for trace, !Dk is the covariant derivative in 
the adjoint representation £ind mo is a mass pcirameter. 

At short distances and small coupling, the derivative term in f(B) 
dominates and Eq.(2.5) reduces to Eq.(2.2). At large distances, the mass 
term dominates and f(B) in Eq.(2.5) reduces to Eq.(2.4) with 

^t = 2ing2 • (2.6) 
Here, we have replaced the parameter mo by a renormsilized pcirameter m. 
In integrating out short-distance degrees of freedom, one expects mo to be 
renormalized. In peirticulcir, m contciins contributions related to the energy 
per unit length of the Wilson lines which enter Eq.(2.5), as we now expleiin. 

An explicit formula for the kernel in Eq.(2.5) for SU(N) is 

-'^ dx 



iV2 



TUX 



X 



j rn^X(x)1 expf- |jda(x'(a) + m^ 



X(0)=y'- 
X(-c)=x 



X 



(2.7) 



P exp i fdaA,x\a) 
x<-y V •'o 



P exp ifdaA,X(a) 

jtiS2 [y^^ ^ ' 



.t2SlJ 



where p denotes the path-order product along X(o) and Ai = X^A^ (X(o)). 
The measure 'dX in Eq.(2.7) is the Feynman one [19] for a particle of unit 
mass in d dimensions. The subscripts si, ti, S2, t2, which run from 1 to N 
for SU(N), are matrix indices: 
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Tr^By(x) (_ ^-^^+mlf\^,y^ By(y)) = 

-1/2 

B\I^H^){-^'''D^+ml) (x,y)B^2t2(^) . (2.8) 

-1/2 

Eq.(2.7) shows that (- ©"^©j^ + nip) (x,y) involves a sum over paths of 

a Wilson line in the adjoint representation which goes from the space-point 
y to the space-point x. Each path contributes with a particular probability. 
As mo is increased, paths of smaller size are weighted comparatively more 
than paths of larger size. The Wilson lines in Eq.(2.7) render Eq.(2.5) gauge 
invariant. When the kinetic term g^jd^x E^ E^(x)/2 in H acts on the Wilson 

line, a contribution to the energy proportional to the length of the non-back- 
tracking p£irt of the path is produced. Hence, it is energetically favorable to 
have mo be non-zero. In the case of a U(l) group, the matrix factor in curly 
brackets in Eq.(2.7) is replaced by 1 and no such Wilson line contribution 
cirlses cind it is energeticcilly favorable for mo to be zero. 

Besides having the correct strong and weak coupling limits, there is 
numericcil evidence from lattice simulations that ^o in Eqs.(2.1) £ind (2.5) is 
a reasonable approximation for the 2+ 1 dimensioneil Yeing-Mills theory. 
Consider expanding f(B) in inverse powers of mo. The result is 

f(B) = ^10 jd 'xTr(B(x)B(x)) + jd 'xTr ( ® ^^Blx) © ^B(x)) + . . . , (2.9) 
where 

Ref.[13] has performed Monte Carlo simulations of the ground state 
function£il at intermediate couplings for SU(2) and found that the form in 
Eq.(2.9) fits the data well with ^io=(0.91 ± 0.02)/g4 and \i2= - (0.19 ± 0.05)/g8. 
Analj^ic strong coupling lattice computations in the Hamilonisin formulation 
lead to simil£ir results.[16] As predicted from Eq.(2.5), ^2 should be negative 
cind this is borne out in simulations. The Monte Ccirlo data implies 
mo» 1.55g2 for SU(2). 
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m. Consequences of The ^proximate Vacuum Functional 

Let us first consider D-2+1 space-time dimensions. In this case, the 
strong-coupling functional in Eq.(2.5) leads to confinement because the 
dimensionally reduced effective theory is similar to a localized Yang-Mills 
theory in two dimensions. The string tension a is obtained from the vacuum 
expectation vEilue of a spaticilly oriented Wilson loop. The result is 

g^mCf 

a= ^ . (3.1) 
The 0"'"''' gluebcill mass Mg is obtained as the coefficient of the 
exponenticil falloff of the correlation function <^By By (0,0) B^^(0,t)> ~ 
c(t) exp(-MGt) for l£irge t. By exploiting Lorentz inv£irl£ince, Mq Ccin be 
extracted from the vacuum functional via 

<^0 1 i By B?j(x) ^ Bji Bji(y) I ^o> ~ c(r) exp(-MG r) , (3.2) 

where r = Ix - yl . In two space dimensions cind at large distances, the 
propagator for B^g in the effective field theory governed by I ^ is 

(-32 + m2)^/2. One finds, using perturbation theory in this effective theory, 

that Eq.(3.2) holds with 

MG»2m . (3.3) 
Eq.(3.3) has a physiccil interpretation. The pciremieter m can be thought of 
as a constituent mass for a gluon. More precisely, m is the effective mass of 
cin adjoint-representation configuration of gluons inside a bound state. Since 
two adjoint representations cire needed to m£ike a singlet, the glueball mass 
is approximately 2m. 

Combining (3.1) and (3.3), one arrives at the main result in Eq.(1.2). 

Monte Ccirlo simulations of the 2+ 1 dimension£il SU(2) Yeing- Mills 
theory have accurately [20] determined the string tension to be 
OMC = (0. 1 12 ± 0.002) g4. Using this value for o in Eq.(1.2), we obtain 
Mg ~ 1.2 g2. The Monte Ccirlo Ccilculations of ref.[21] give Mg 
= (1.59 ± 0.01) g2. Thus, we find MQ/'\j c =3.6, while numerical simulations 
[21] give Mg/V^= 4.77 ± 0.05. Using the difference between our results and 
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those of the lattice as a mecins of estimating systematic uncertciinty, our 
method for computing the O"*""*" gluebcill mass is accurate to about 25%. For a 
summEiry of computations of Mg/V^ > see the references in ref.[22]. 

For SU(3), Monte Carlo simulations give omc = (0.307 ± 0.004) g4.[23] 
This value of o in conjunction with Eq.(1.2) leads to the SU(3) results 

MG»1.84g2 , Mg/Vo' = 3.3 . (3.4) 
The result in Eq.(3.4) is in reasonably agreement with strong coupling 
Hamiltonian computations [24]. 

Recently, another analytic approach, based on finding exact 
eigenstates of the kinetic energy term in H, gives Mq - N/n g2 for SU(N).[25] 
For SU(2), one obtains Mg = 0.64 g2, which is considerably smaller than the 
Monte Carlo result, suggesting that corrections involving the potential 

energy term jd^x B^2(^^/^S^ may be importcint. For SU(3), the approach 

gives Mg ~ 0.95 g2, which is about heilf the veilue in Eq.(3.4). 
IV. Extrapolation to D=4 

The result in Eq.(1.2) does not directly apply to D=3+l dimensions. 
However, one can appeal to the idea of dimensional reduction. [26] If a 3+1 
dimensional gauge theory is confining, then the computation of spaticil 
correlation functions cind Wilson loops involves only those degrees of 
freedom within the order of Lc in the time direction, where Lc is the 
confinement length. Hence, it should be possible to approximate a 3+1 
dimensional confining theory by a 2+1 dimension£il theory. This idea has 
been used several times in past work[ 15,27] £ind is embodied in the result 
that fsc(B) Ccin be used for the purposes of computing l£irge spatial Wilson 
loops. Due to the beta function being negative, |j. is expected to be 
exponentially large in l/g2. In short, assuming four-dimensional 
confinement, results for Mg/V^ should be approximately the same in 2+1 
and 3+1 dimensions. In fact, Monte Carlo simulations of the SU(2) gauge 
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theoiy indicate that the glueball spectra, in units of y a, in 2+ 1 £ind 3+ 1 
dimensions agree to about 15%. [21] 

Assuming the v£ilidity of dimensioneil reduction, Eq.(3.4) with 
^[c = 440 MeV gives a v£ilue of the 0"^"^ gluebcill mass Mq in D=3+ 1 
dimensions for the SU(3) gauge theory of 

MG«1.5GeV . (4.1) 

V. Summary 

In summEuy, by using £in approximate vacuum function£il which 
interpolates between strong and weak coupling forms, we have obtained the 
relation in Eq.(1.2) between the lightest O"*""*" glueball mass and the string 
tension in 2+ 1 dimensions. Using dimensional reduction, we obtciin a v£ilue 
of about 1.5 GeV for Mg in 3+1 dimensions. 

Our value of 1.5 GeV for Mq is much larger than many of the natural 
scciles in Yang-Mills theory: ^fa is 440 MeV, the deconfining phase transition 
temperature is -250 MeV, the mass of the color-singlet gluon cloud around 
a qucirk is -300 MeV assuming that the contribution to the constituent mass 
of a quark in a bound state comes from such a gluon cloud, and the 
topological susceptibility <v^> which enters the r\' mass is <v2>i/4 » 180 MeV 

[28-33]. On this basis one might expect approximate anal5^ic computations 
of Mg to yield values less than 1 GeV. For example, if the continuum strong 
coupling 2+1 dimensional Hamiltonian result [25] is assumed to extrapolate 
to D=3+ 1 by dtmensioneil reduction, one would obtain Mg ~ 760 MeV. 

Since, as mentioned above, the error in our results in 2+1 dimensions 
is estimated to be about 25% cind that the error in extrapolating to 3+1 
dimensions is about 1 5%, the totcil error in Mg is around 30%. Our VEilue of 
1.5 GeV for Mg is about 15% less thsin the most accurate current lattice 
value of 1740 MeV [6]. 
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